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Abstract 

We use the gauge-gravity duahty conjecture to compute spectral functions of the stress-energy tensor in 
finite temperature A/'=4 supersymmetric Yang-Mills theory in the limit of large Nc and large 't Hooft 
coupling. The spectral functions exhibit peaks characteristic of hydrodynamic modes at small frequency, 
and oscillations at intermediate frequency. The non-perturbative spectral functions differ qualitatively from 
those obtained in perturbation theory. The results may prove useful for lattice studies of transport processes 
in thermal gauge theories. 
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I. INTRODUCTION 



Theoretical analysis of the properties of strongly interacting hot and dense matter is a hard 
problem. Even when the density of baryons in thermal equilibrium is negligible, perturbative QCD 
calculations are only reliable at temperatures which are much higher than the temperature of 
the deconfinement transition P|. While lattice simulations can provide an account of equilibrium 
thermodynamic properties of the theory, questions involving real-time dynamics such as momentum 
transport, thermalization, and various production rates are much harder to answer. For near- 
equilibrium states, physical information is contained in equilibrium response functions: for example, 
the dilepton rate is proportional to the spectral function of vector currents, and the viscosities are 
determined by the spectral function of the relevant components of the stress-energy tensor. 

Thus it is valuable to study models where analytic results for the real-time response functions 
can be derived. A number of recent studies focused on a particular model, the AA=4 supersymmet- 
ric SU{Nc) Yang-Mills (SYM) theory at finite temperature. The interest in this theory is due to 
Maldacena's gauge-string duality conjecture 0] which provides an effective description of the the- 
ory's non-perturbative regime in terms of semiclassical gravity in a five-dimensional asymptotically 
Anti de Sitter space. 

The SYM theory is conformal, and has only one tunable parameter, the 't Hooft coupling A. 
Thermal equilibrium is characterized by the black-body equation of state at any non-zero tempera- 
ture, and can not (at large 't Hooft coupling) be viewed as a gas of weakly interacting quasiparticles. 
The standard hydrodynamic singularities in two-point functions of conserved currents in SYM the- 
ory were found from the dual gravity description in the limit of large A and large Nc Q, |^ . The ratio 
of the shear viscosity to volume entropy density in this theory has been found to be rj/s = l/47r 

Q], which is a much smaller number than the corresponding result in a weakly coupled gauge 
theory. 

In this paper, we focus on the full spectral functions (rather than their small-frequency limit) of 
the stress-energy tensor in the SYM theory at large 't Hooft coupling and large Nc- The knowledge 
of the full spectral function is important for two reasons. On the one hand, our results provide 
the first example of a non-perturbative spectral function calculation in a strongly coupled four- 
dimensional gauge theory at finite temperature, obtained without performing lattice simulations. 
On the other hand, the non-perturbative spectral functions of SYM theory may be useful for lattice 
computations of transport coefficients in realistic gauge theories at temperatures not too far from 
the deconfinement transition. Indeed, a Euclidean correlation function (computed on the lattice) is 
proportional to the integral of the real-time spectral function over all frequencies. In the simplest 
method of reconstructing the spectral function from the Euclidean data, one first assumes an 
ansatz for the spectral function, and then fits the parameters of the ansatz to the lattice data 
This method has been applied to the computation of the shear viscosity in QCD 8]. Theoretical 
understanding of what the correct non-perturbative ansatz might be is of primary importance for 
this approach. Our results for the non-perturbative spectral functions of SYM theory may therefore 
prove useful for lattice studies of transport processes in thermal gauge theories. 
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II. CORRELATION FUNCTIONS FROM GAUGE-GRAVITY DUALITY 



The spectral function Xfiv,ai3{k) is defined as 

X>.uMk) = f d'x e-'''mT^,{x),T^pm) . (1) 



It is proportional to the imaginary part of the retarded Green's function, XiJ.u,ai3{k) 
—2linG^u,ai3{k) , where 



G^,,ap{k) = -iJ d''xe-'"'^0{x^){[T^,{x), r„^(0)]) , (2) 

and is an odd, real function of k. The shear and bulk viscosities are proportional to the zero- 
frequency slope of the specific components of the spectral function, for example, 

V = ^^^^Xxy,xy{k'^,k=0) . (3) 

The retarded correlation function of the stress-energy tensor of SYM theory admits a simple decom- 
position (we follow the notational conventions of Ref. Q). At zero temperature, Lorentz symmetry 
combined with conservation and tracelessness of T^i/ implies that has the form 

where H^^^ap = \ {P/MaPup + PfilsPua) ~ '^P^iuPaf) is s projector onto conserved traceless symmetric 
tensors, P^u = i^^y — k^ky/k"^, and k"^ = —{k^)"^ + k"^. At non-zero temperature, G^i,^ai3 in a 
conformal theory can be described by three symmetry channels 

Gf^u,al3{k) = Sfj_u,a(3 Gi + Q^iy,al3 G2 + L^^^a/i G^ , (4) 

where S^j^u^af}-, QfMu,ai3, Lf^^^a/s are the appropriate orthogonal projectors which provide three inde- 
pendent Lorentz index structures 0]. Choosing the spatial momentum along x^, k^ = {—uj, 0, 0, q), 
the components of the correlation function are 

1 g2 

Gtx\tx^{k) = — oGii^^^q), (5a) 
2 

Gtt,tt{k) = — 2^^2(0;, , (5b) 
Gx^x^,x^x^{k) = ]^G^{uj,q) , (5c) 

with all other related to the above by the rotation invariance. As a function of complex w, in the 
low-frequency limit, Gi{uj,q) has a shear-mode singularity, G2{uJ,q) has a sound-mode singularity, 
and G3{u;,q) has no hydrodynamic singularities. In the limit of vanishing 3-momentum, Gi{uj) = 
G2{ijj) = G^iuj). At zero temperature, Gi = G2 = G3 = Gs- 

In the regime of large 't Hooft couplin g, t he three scalar functions Ga{ijJ-,q) can be computed 
using the gauge-gravity duality recipe jfll. Ilfl|. The duality essentially reduces the computation 
of a two-point correlation function to a boundary- value problem for a linear ordinary differential 
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equation. For the zero-temperature theory, the retarded two-point function can be found for 
example in 0]^, 

Gs{k) = (ln|fc2| _^^^(-fc')sign^) . (6) 

In order to compute the retarded correlators at non-zero temperature, one has to analyze the 
"wave equations" (one for each symmetry channel) which describe propagation of the corresponding 
metric perturbations in the AdS-Schwarzschild background spacetime of the dual description. The 
differential equations are of the form 

d 

-^Za{u) +Pa(u) — Za{u) + qa{u)Za{u) =0, (7) 

where the coefficients Pa{u), qa{u) [to be specified shortly] depend on the dimensionless frequency 
rx> = u)/2'kT and momentum q = g/27rT, and a = 1,2,3 labels the three symmetry channels. The 
coordinate u ranges from to 1, where u = corresponds to the boundary of the asymptotically 
AdS spacetime, and n = 1 corresponds to the event horizon of the background. 

For all three equations (jT)), the characteristic exponents at u = are equal to and 2, and the 
exponents at u = 1 are ibii»/2. Information about the retarded correlation function is encoded in 
the solutions to Eq. Q which satisfy the incoming wave condition at the horizon, corresponding to 
the exponent —ivo/2 at u = 1. The correct solution is thus of the form Za{u) = (1 — u)^'^^^'^Fa{u), 
where Fa{u) is a regular function at the horizon. The solution satisfying the incoming- wave con- 
dition at the horizon can be written as a linear combination of two independent local solutions at 
li = 0, 

Za{u) = AaZiiu) + BaZi\u) , (8) 

where Zl{u) and Z^/ {u) are given by their standard Frobenius expansions as 

Zi = l+ fe^y u + ha Zi\u) In n + h^^} + ■ ■ ■ , (9) 
Zi^ = u\l + b^lu + b^Sin' + ■■■). (10) 

All the coefficients h^^} jj (except h^^)) and ha are determined by the recursion relations obtained 
by substituting the above expansion in the differential equation ((T)). The coefficient is left as 
a free parameter, reflecting the fact that one can always redefine the local solutions by adding a 
constant multiple of Zl^{u) to Za{u). Without loss of generality, we set b^^j = thus fixing the 
definition of Za^\u). The retarded functions Ga are then given by^ 

G„(c.,g) = -vr2iV,2r^MM). (11) 

•^a\^, Q) 



^ At zero temperature, the dual gravity result for the correlator in strongly coupled A/'=4 SYM theory coincides 
with the one obtained in free field theory, due to a non-renormalization theorem IlLl . 

The real parts of the correlator in the shear and sound channels have additional contributions arising from writing 
the gravitational action in terms of the gauge-invariant variables 0. These terms are irrelevant for the present 
discussion. 
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FIG. 1: Finite-temperature part of the spectral function for transverse stress {xxy,xy — xiyxy)/^: plotted 
in units of tt'^N^T^ as a function of dimensionless frequency to = Lo/2nT. Different curves correspond to 
values of the dimensionless spatial momentum q = q/2TTT equal to 0, 0.6, 1.0, and 1.5. 



As is evident from Eqs. ©, (|1U|) . the coefficient A is given by the boundary value of the solution, 
Aa = lim^^o Za{u), while the coefficient B can be expressed in terms of the boundary value of the 
second derivative of the solution, 



Ba = l lim - 2AahaHu)) - ^Aaha • 



(12) 



From the recursion relations one determines ha = —5(9^ — ro^)^, which is an analytic function of 
K), q, and therefore represents a contact term which we drop. The retarded correlator is therefore 
equal to 



Ga = -vr^iV^T^ lim ( _ ha ln(u] 



(13) 



III. SPECTRAL FUNCTIONS OF THE STRESS-ENERGY TENSOR 

In order to determine the spectral function for transverse stress, we need to solve the master 
equation Q whose coefficients are given by 0] 



1 + 
uf 



up 



(14) 



where f = \ — v? . Solving Eq. numerically, we find Xxy,xy{^, l) as explained in the previous 
section. The spectral function is linear in to for small frequencies, Xxy,xy = 
0(t»2,q2)), then increases monotonically, and at large frequencies it asymptotes to the zero- 
temperature result Xxy',xy — f — to^)^ 7r^N^T^9{—k'^). When this zero-temperature contribution 
is subtracted, the resulting function exhibits oscillations which damp rapidly as frequency grows. 
The oscillations appear around u) = q, and their amplitude grows with q. The unsubtracted Xxy,xy 
is positive, as it should be. FigureHshows graphs of Xxy,xy for several values of three- momentum^ . 



The numerical procedure outlined above allows one to compute both real and imaginary parts of Gpi^^Q/S • The real 
part which can in principle be reconstructed from the spectral function, also exhibits oscillations which damp as 
frequency grows 



n p i 
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FIG. 2: Left: spectral function for longitudinal momentum density, Xtx,tx, plotted in units of tt^A^^T^, as a 
function of dimensionless frequency to = lojlixT . Different curves correspond to values of the dimensionless 
spatial momentum q = q/2'KT equal to 0.3, 0.6, 1.0, and 1.5. At large to, the curves asymptote to the 
zero-temperature result ^q^(ro^ — q^). Right: spectral function for energy density, Xtt.tt, plotted in units of 
tt'^N^T'^, as a function of dimensionless frequency vo = uj/2TrT. Different curves correspond to values of the 
dimensionless spatial momentum q = g/27rT equal to 0.3, 0.6, 1.0, and 1.5. At large i», the curves asymptote 
to the zero-temperature result ■|4q''/3. 



The shear viscosity follows from the Kubo formula r] = ttN^T^/S, and is in agreement with 
the earlier results Q| . 



In the shear channel, the coefficients of the master equation are given by 

The spectral function Xtx,tx{^T q) is shown in Fig. |21 At small momentum, the spectral function 
exhibits a narrow peak at to = q^/2, characteristic of the hydrodynamic shear mode. 
In the sound channel, the coefficients of the master equation ((Tj) are given by 0| 

M/(3n)^ -I- q^(n^ — 3)) 

, , 3«)4+q4(3_4^2^^4) ^q2(4^2^2_g^2_4^3j) 

uj^{3x»^ + q''{w^ — 3)) 

The spectral function xtt,tt{^,Q) is shown in Fig. |2 At small momenta, the spectral function 
exhibits a narrow peak at to = q/\/3, characteristic of the hydrodynamic sound mode. For both 
Xtx,tx and Xtt,tt, the finite-temperature contributions have oscillatory behavior, similar to the one 
seen in Xxy,xy 



IV. DISCUSSION 



It is easy to understand the limiting behavior of the spectral functions found above. The lo, 
q dependence at small frequency is predicted by the linearized hydrodynamics, while the large co 
dependence is fixed by the scale invariance of the SYM theory. An intriguing feature is the presence 
of oscillations in the finite-temperature contribution, which appear around uj = q, and then decay 
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rapidly. Mathematically, such damped oscillations are due to the characteristic asymptotic behav- 
ior, ~ exp (—am), of the solutions to Eq. (O, where a is a complex number 0, Q|. They reflect 



the presence of an infinite sequence of poles in the lower half-plane of the retarded correlators 

One can compare the spectral functions of the strongly coupled SYM theory with the pertur- 
bative results in a weakly coupled scalar or pure gauge theory discussed in Ref. Q]. At weak 
coupling, the spectral function for transverse stress at ^=0 grows linearly at small frequency, and 
is proportional to at asymptotically high frequency. In between, however, there is a range where 
the spectral function decreases as 1/uj. This behavior at intermediate frequencies at weak coupling 
differs from what one observes in the strongly coupled SYM theory, where the spectral function 
grows monotonically for all uj. Based on the numerical results for X^j.u,ai3{^jQ): one may propose 
the following ansatz for the zero-momemtum spectral function for transverse stress 

Xxy,xy{uj) = Xly,xyi^) + Im ^ Q e'^^ ^ , (18) 

i 

where q, are real coefficients, and Re(aj) > (assuming positive frequency). 

In this paper, we computed the spectral functions in the simplest thermal gauge theory with a 
known gravity dual. It should also be feasible to find spectral functions for a class of non-conformal 
gauge theories, such as the one analyzed in Ref. 

Note added: while this paper was being completed, we became aware of the work on the same 



subject by D. Teaney |l7|. 
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